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We report the first observation of the radiative decay B° — > (j)K°Y using a data sample of 
772 x 10 6 BB pairs collected at the T(45) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. We observe a signal of 35 ± 8 events with a significance of 
5.4 standard deviations including systematic uncertainties. The measured branching fraction is 
.3§(B° <j)K°y) = (2.66 ±0.60 ±0.32) x 10~ 6 , where the uncertainties are statistical and sys- 
tematic, respectively. We also report the first measurement of time-dependent CP violation pa- 
rameters: J^o r = +0.74+? : g(stat)±H : '°(syst) and = +0.35 ± 0.58 (stat)+°;^(syst). We 
also precisely measure 3§{B + -> <$>K + y) = (2.34±0.29±0.23) x 10~ 6 . The observed M$ K mass 
spectrum differs significantly from that expected in a three-body phase-space decay. These results 
are preliminary. 
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Rare radiative B meson decays play an important role in the search for physics beyond the 
standard model (SM). These decays are forbidden at tree level in the SM, but allowed through elec- 
troweak loop processes. The loop can be mediated by non-SM particles, and therefore is sensitive 
to new physics (NP). Exclusive b — > sy decays have been extensively measured, but their total sum 
so far accounts only for 44% of the inclusive rate. Therefore, further measurements of branching 
fractions for exclusive B — > <pKy modes will improve our understanding of the b — > sy process. The 
emitted photons are predominantly left-handed (right-handed) in b — > sy (b — > sy) decays. This 
suppresses the CP asymmetry in the SM by the quark mass ratio (2m s /mb). The expected mixing- 
induced CP asymmetry parameter {5^) is ^(3%) and the direct CP asymmetry parameter (s>/) is 
~ 0.6% [[]]]. However, in several extensions of the SM, both the photon helicities might contribute 
to the decay. Therefore, any significantly larger CP asymmetry would be clear evidence for NP. 

In this presentation, we report the first observation of neutral mode B° — > (j)K°y the first 
measurements of time-dependent CP violation in this mode, as well as improved measurements in 
B + — > <pK + y using a data sample of 772 x 10 6 BB pairs collected at the T(45') resonance with the 
Belle detector ^ at the KEKB asymmetric-energy e + e~ (3.5 on 8.0 GeV) collider This data 
sample is nearly eight times larger than the sample used in our previous measurement [§J. In the 
decay chain T(45) — > B°B° — > f rec J \ ag , where one of the B mesons decays at time ? rec to the signal 
mode / rec and the other decays at time ? tag to a final state / tag that distinguishes between B° and B°, 
the decay rate has a time dependence given by 



&(At) = <! 1 +q y sin ( Am d At) + £</ cos (Am d At) 

4t b o 



(1) 



Here Tgo is the neutral B lifetime, Am d is the mass difference between the two neutral B mass 
eigenstates, At = ? rec — h ag , and the ^-flavor charge q equals +1 (—1) when the tagging B meson is 
a B° (B°). 

Signal candidates are reconstructed in the B + — > <j)K + y and B° — > (pK^y modes, with (j) — >■ 
K + K and — > n + n~. The invariant mass of the kaon pairs from is required to be within 
\Mk+k- — 1 < 0.01 GeV/c 2 , where m§ denotes the (j) meson world-average mass The K® 
selection criteria are the same as those described in the Ref . [|7|] ; the invariant mass of the pion pairs 
is required to satisfy M n + n - G (0.482,0.514) GeV/c 2 . The high energy prompt photons are selected 
from isolated clusters within the barrel region of the calorimeter and center-of-mass system (cms) 
energy, Ey ms 6 (1.4,3.4) GeV. B candidates are identified using two kinematic variables: the energy 
difference AE = Eg ms — E^ m and the beam-energy-constrained mass M\, c = \/ (-^beam) 2 — (^s ms ) 2 ' 
where E^ m is the beam energy in the cms, and £^ ms and /?| ms are the cms energy and momentum, 
respectively, of the reconstructed B candidate. In the M\, c calculation, the photon momentum is 
replaced by (E^ m — E^™) to improve its resolution. The candidates that satisfy the requirements 
Mbc > 5.2 GeV/c 2 and |A2s | < 0.3 GeV are selected for further analysis. 

The dominant background is from the continuum process, which is suppressed by a require- 
ment on likelihood ratio {&) from event shape variables and the B flight direction. We require 
Si > 0.65, which removes 91% of the continuum while retaining 76% of the signal. In the 
B° -> (j>K$Y mode, some backgrounds from b — > c decays, such as D°7i°, D 0r q and D p + peak 
in the M\, c distribution. We remove the dominant peaking backgrounds by vetoing <pK® combina- 
tions consistent with the nominal D mass [Eft. Some of the charmless backgrounds, where the B 
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Figure 1: The AE and M bc projections for B + ->• <j)K+y ((a) and (b)) and B° ->■ i^y ((c) and (d)). The 
points with error bars represent the data. The curves show the total fit function (solid red), total background 
function (long-dashed black), continuum component (dotted blue), the b — s> c component (dashed-dotted 
green) and the non-resonant component as well as other charmless backgrounds (filled magenta histogram). 



meson decays to 0^(892), (j)KK° and (j)Kri also peak in Mb c , but are shifted towards lower AE. 
Another significant background is non-resonant B —¥ K + K Ky, which peaks in the AE-Mbc signal 
region; it is estimated using the mass sideband M K + K - G (1.05, 1.3) GeV/c 2 , in data. 

The signal yield is obtained from an extended unbinned maximum-likelihood (UML) fit to the 
two-dimensional AE-M\, C distribution as shown in Fig. [| The probability density functions (PDFs) 
for different components and fitting procedure are described in detail elsewhere The fit yields 
a signal of 136 ± 17 B + — > §K + y and 35 ± 8 B° — > tyK^y candidates. The signal in the charged 
mode has a significance of 9.6a, whereas that for the neutral mode is 5.4a, including system- 
atic uncertainties. The observed (j)K mass spectrum differs significantly from that expected in a 
three-body phase-space decay. Nearly 72% of the signal events are concentrated in the low-mass 
region {Mqk G (1.5,2.0) GeV/c 2 ). The Monte Carlo (MC) efficiency after reweighting accord- 
ing to this Mtj,K dependence is (15.3 ±0.1)% for the charged and (10.0 ±0.1)% for the neutral 
mode. We measure the branching fractions as 38{B + -s- (j)K + y) = (2.34 ±0.29 ±0.23) x 10~ 6 and 
38(B° (j)K°y) = (2.66 ±0.60 ±0.32) x 10~ 6 , where the uncertainties are statistical and system- 
atic, respectively. 

For the CP asymmetry fit, we select events that satisfy M\, c G (5.27,5.29) GeV/c 2 and AE G 
(—0.2,0.1) GeV. The vertex position for the / rec decay is reconstructed using the two kaon tracks 
from the (f) meson and that of / tag decay is obtained from well-reconstructed tracks that are not 
assigned to f mc We then use a flavor tagging algorithm [ |l0| ] to obtain q and the tagging quality 
factor r G (0,1). The typical vertex reconstruction efficiency (z resolution) is 96% (115/xm) for 
/ rec and 94% (104 /im) for f tag . After all selection criteria are applied, we obtain 75 (436) events 
for the CP fit with a purity of 45% (35%) in the neutral (charged) mode. The signal PDF is 
given by a modified form of Eq. ([]]) after fixing Tgo and Amj to their world-average values ^ and 
incorporating the effect of incorrect flavor assignment. Since the non-resonant component has the 
same NP as the signal B — > (j)Ky, we treat this as signal. The only free parameters in the CP fit are 



y and srf, which are determined to be = ±0.74+ t 1 ' ^(stat)lJ] ^(syst) and 



±0.35 ± 



0.58(stat)_Q j () (syst). We define the raw asymmetry in each At bin by (/V+ —N-)/(N + +N-), 
where (A^_) is the number of observed candidates with q = ±1 (—1). Figure || shows the 
At distributions and raw asymmetry for events with good tagging quality (r > 0.5, 48% of total). 
We find that the error on y in the MINUIT minimization is much smaller than the expectation 



from MC simulation and has a probability of only 0.6% [11]. This is due to low statistics and the 
presence of a single special event. Instead, we use the ±68% confidence intervals in the residual 
distributions of 5? and srf, determined from the toy MC as the uncertainties for the final result. 
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Figure 2: Af distributions for q = + 1 and — 1 (left) and the raw asymmetry (right) for well-tagged events. 
The dashed curves in the Af plot are the sum of backgrounds while the solid curves are the sum of signal and 
backgrounds. The solid curve in asymmetry plot shows the result of the UML fit. 

In summary, we report the first observation of radiative B Q — > <pK°y decays using a data sam- 
ple of 772 x 10 6 BB pairs. The observed signal yield is 35 ± 8 with a significance of 5.4a. We 
also report the first measurements (preliminary) of time-dependent CP violation parameters in this 
mode: y = 0.74±°jg+ n i° and srf = 0.35 ±0.58+^. We precisely measure the branching fraction 
and charge asymmetry (£/cp = —0.03 ± 0. 1 1 ± 0.08) for the charged mode. The signal events are 
mostly concentrated at low <pK mass, which is similar to a two-body radiative decay. With the 
present statistics, these measurements are consistent with the SM predictions and there is no indi- 
cation of NP from right-handed currents in radiative B decays. Much more luminosity is necessary 
for a precise test of the SM. 

The author wish to thank the KEKB group for excellent operation of the accelerator, the KEK 
cryogenics group for efficient solenoid operations, and the KEK computer group and the Nil for 
valuable computing and SINET3 network support. 
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